Preprint typeset in JHEP style - HYPER VERSION 



Holographic Spectral Functions in Metallic 
AdS/CFT 



Javier Mas? Jonathan P. Shock and Javier Tamo* 

Departamento de Fisica de Particulas, Universidade de Santiago de Compostela 
and 

Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E-15782 Santiago de Compostela, Spain 



Abstract: We study the holographic D3/D7 setup dual to M = 4 supersymmetric 
Yang-Mills with quenched fundamental matter. We extend the previous analyses of 
conductivity and photoproduction to the case where there is a finite electric field. Due 
to the electric field a special region in the D7-brane geometry, labelled the singular 
shell, appears generically, and the computation of correlators involves a careful study 
of the indicial exponents both at this singular region and at the horizon. We show 
that there is a unique choice consistent with the known expression for the electrical 
conductivity found by Karch and O'Bannon [1]. We explore the parameter space 
spanned by the quark mass, the baryon density and the electric field. We find a 
region where the conductivity and photoproduction change rapidly and trace this 
behavior to competing effects which manifest themselves as a crossover behavior in 
the probe brane embeddings. 
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1. Introduction 

Since its original formulation, the AdS/CFT correspondence has been applied to 
an extremely diverse range of phenomena. Only recently, the applications to non- 
relativistic systems [2] and superconductivity [3, 4] have been elucidated and a large 
number of papers on these subjects continue to be published. One of the most recent 
applications of the correspondence has been to study gauge theories in the presence 
of electric and magnetic fields [1, 5-12], a useful way to explore the properties of 
many physical systems. 

In this paper we are interested in the study of strongly coupled plasmas, in partic- 
ular the Dp/Dq brane intersections holographically dual to Yang-Mills theory with 
fundamental hypermultiplets [14, 15]. We will use the benchmark D3/D7 system in 
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order to study the properties of such plasmas, generalising wherever possible to the 
Dp/Dq intersection. 

One of the first additions to the most basic AdS/QGP (Anti-de-Sitter/Quark Gluon 
Plasma) calculations was the introduction of a finite charge density [16, 17]. This 
corresponds to studying a Dg-brane probe in a black Dp-brane background (a geom- 
etry with an AdS boundary and a black hole at the centre) with the time component 
of a world-volume gauge field on the .Dg-brane turned on. One of the most striking 
consequences of the addition of finite charge density was the disappearance of the, so 
called, "Minkowski embeddings" [17]. These Minkowski embeddings correspond to 
.Dg-branes which wrap vanishing cycles which end away from the black hole horizon. 
These solutions are characterised by a discrete spectrum of excitations corresponding 
to infinitely long-lived bound states. However, no matter how small the quark den- 
sity introduced, the bulk gauge field flux lines need somewhere to end, and the only 
option is the horizon, removing the possibility of having Minkowski embeddings. An 
immediate consequence of this is that the introduction of heavy flavors with finite 
charge density leads to flavor probe brane embeddings which, as they run in from 
the AdS boundary into the IR, are flat all the way down to some small region where 
they develop a long spike which eventually hits the horizon. All embeddings in this 
setting then take the form of "Black Hole" , or "Schwarzschild embeddings" . In more 
physical terms, this means that with any finite quark density, all mesons are unstable 
and the spectral function becomes continuous. It also means that the "fundamental 
phase transition" between stable and unstable mesons disappears. 

In this setting, the magnitude of the charge density can be used as a control pa- 
rameter in which the smooth, continuous spectral functions tend, in the vanishing 
charge density limit, to those of the stable phase with its signature discrete spectral 
function. In this way one can approach the stable meson phase arbitrarily close. In 
the limit where the spectral function takes the form of sharp peaks one can study 
the dispersion relation of the modes in the plasma by tracking the peak positions in 
the frequency plane as a function of their momenta. This was performed for trans- 
verse excitations in [18] and for the more involved longitudinal ones in [19]. For the 
transverse case subluminal asymptotic velocities were observed, while in the second, 
some evidence for superluminal propagation was obtained. This second, rather more 
exotic behavior is currently under study by looking at the fully fledged quasinormal 
mode calculation [20]. 

The next step in exploring these kinds of systems is the introduction of electric 
and/or magnetic fields within the holographic setups at zero and finite temperature. 
In [1, 7-9, 21] the phase structure of these systems was studied in detail, while 
the spectrum of excitations was also looked at in a certain subsector of the phase 
diagram. In the present study we will be interested in extending this work to the 
rest of the phase diagram in the case of the electric field, where the behavior is more 
complicated than in the magnetic case. For similar studies in the Sakai Sugimoto 
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model see [11-13]. 

In the electric case, for small enough ratios of the quark mass to the electric field 
strength, the probe Dg-brane becomes unstable, the signature of this instability 
being a region on the brane's worldvolume where the lagrangian becomes complex 
[1, 9]. The solution to this pathology is the introduction of a current in the direction 
of the electric field. Though this removes the pathology, the remnant of this is a 
submanifold of the Dg-brane worldvolume where the lagrangian density vanishes on- 
shell. This submanifold is known in the literature variously as the "singular shell" 
[8] and the "ergosphere" [10] (due to its interpretation in the T-dualised geometry). 
It turns out that for any given electric field strength, there is a unique current which 
will cure the pathology. The ratio of the current to the electric field then provides 
a definition of the D.C. conductivity, as would be computed via a simple Ohm's law 
experiment. 

Although the singular shell cures the problem of having an ill-defined action, it does 
introduce new complications. In particular in order to compute the spectral function 
for excitations on top of classical embeddings which pass through the singular shell, 
the boundary conditions for the excitation, which are non-trivial, must be under- 
stood. In this paper we solve this problem, thereby allowing us to calculate more 
physical quantities in this interesting setup. 

A check that we have the correct boundary conditions will be to compare the D.C. 
conductivity calculated from the macroscopic setup [1] to that calculated using the 
Kubo relation, which involves the excitations on top of the classical brane embedding. 
Once the correct boundary conditions have been found we will be able to calculate 
the spectral function for light-like momenta and thereby study the photoproduction 
rate for the plasma in the presence of an electric field. 

The parameter space that we are interested in is spanned by the magnitude of the 
electric field, the quark mass and the charge density, and in the course of studying 
the spectral function we will discover a region of particular interest in it. In this 
region there is a very rapid change in behavior of the spectral function, going quickly 
from smooth oscillations to sharp spikes as we traverse this region of the tunable 
parameters. The change in behavior can be seen in several distinct physical quan- 
tities, including the conductivity, the susceptibility and the photoproduction rate. 
The behavior can be traced back to a rapid crossover in the embedding solutions 
which affects all subsequent quantities. 

The outline of the paper is as follows: In section ||] we introduce our holographic 
setup. Wherever possible we work within the framework of the general Dp/Dq 
flavor brane intersection. In this context we derive the constants of motion and the 
conductivity as given in [1] . 

In section |] we turn to the specific case of the D3/D7 intersection and illustrate 
the various types of brane behavior present by a series of plots. We then examine 
the change in behavior of the profiles upon a continuous variation of the electric 
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field and the quark mass, as described above, and find the region where there is a 
quick crossover in behavior. We also study the effect of this transition on the chiral 
condensate, charge susceptibility, conductivity and diffusion which can be derived 
directly from the background embedding solutions. 

In section |] we move on to the computation of correlators. We do this for transverse 
fluctuations only, propagating in the direction of the background electric field. We 
show that this is enough to compute the conductivity using the Kubo relation. This 
can be done analytically for general Dp/Dq systems and we find perfect agreement 
with the expression in [1]. 

In section [5], in order to calculate the correlation functions away from the hydro- 
dynamic limit we again specialize to the case of the D3/D7 setup. We study the 
transverse correlator at zero momentum for different values of the electric field and 
see that the crossover in the embedding behavior has a large effect on the spectral 
function. From the same function evaluated on the light-cone we obtain the photo- 
production rate. On integrating this magnitude to obtain the total luminosity we 
find a dramatic increase in the total number of emitted photons upon traversing the 
crossover region. 

We finally summarize our results and comment on possible extensions to the current 
work. 



2. Holographic setup: background solutions 

In this section we summarize the framework used to study the embedding of Nf probe 
.Dg-branes in the background of a stack of N c ^> Nf coincident Dp-branes. This has 
been studied in many contexts in the literature, and a review of such investigations 
can be found in [22]. We will start here by detailing the generic Dp/Dq intersection 
and later specialise to the D3/D7 system. 

On a stack of coincident Dq-probe branes there is a natural U(Nf) global symmetry 
whose abelian center can be identified with a baryonic U(l) symmetry. The dynamics 
of the brane embeddings and of the U(l) gauge field is dictated by the probe Dq- 
brane DBI action 

Sdbi = ~N f T Dq [ dx 9+1 e~ V- det(<U + 2na'F ab ), (2.1) 

JD q 

where T Dq = l/((2n\/a') q g s \/a') is the Dg-brane tension, g s the string coupling, 
a' the inverse string tension and g a b the pullback metric from the 10-dimensional 
background. Here we will not deal with isospin degrees of freedom, hence any possible 
non-abelian effects are ignored. The background in which we wish to place these 
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probe branes is the near-horizon limit of a stack of non-extremal Dp-hranes 



dsj = G AB dX A dX B = if-s (-f(r)dx 2 + dx 2 p ) + H* ( + r 2 d{l 2 8 „ p ) ,(2.2a) 



J(r) 

e 4> = #(3-p)/4 ^ ^ = #-1 ^ ( 2>2b ) 



#(r) 



f(r) = 1 - (^) , (2.2c) 
and the Hawking temperature is given by 



T = 7 ~ P 



The probe branes wrap an n-sphere in the directions transverse to the Dp-branes, 
so we write 

dfi 8 _p = d# 2 + sin 2 6dQ 2 n + cos 2 6dQ 2 ^ p _ n , (2.4) 

setting the classical Dg-brane embedding by a functional dependence 1 cos 6 = ip(r). 
The Dg-brane pullback metric is written as 

ds 2 = g ab dx a dx b = g 00 (r)dx 2 + ga(r)dx 2 p + g rr (r)dr 2 + g nn (r)dtt 2 n , (2.5) 

where the components of the induced metric g a b will coincide with the components 
of the background 10-dimensional metric G a b except for the radial term, given by 

g„. = G rr + Geem\rf = (l + r 2 f{r)^^j . (2.6) 

All of the following analytic results will hold for any AdS black hole metric whose 
time component goes to zero linearly at the horizon, independent of the coordinate 
system chosen. All physical results will also be independent of redefinitions of the 
radial coordinate. Here we will consider both finite charge density and a finite elec- 
tric field directed along one of the spatial directions of the Minkowski space, say 
x v . The holographic ansatz for the world-volume gauge field that takes these two 
contributions into account is 



A = A (r)dx° + (E p x° + A p (r)) dx p . (2.7) 

A (r) and A p (r) will be normalizable solutions to the equations of motion, and their 
leading terms in an expansion around the asymptotic AdS boundary correspond to 
the vacuum expectation values (vevs) of current densities J°(x) and J p (x) respec- 
tively. The non- normalizable piece, E p x°, is dual to the source of the current, namely 
the electric field itself. 

1 this is possible by exploiting an 0(8 — p — n) symmetry. 
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For compactness we will define a new matrix given by ^ a b = g a b + {2^ ex!) F a b. The 
only non-diagonal components of this matrix, 7^, are found when i,j take values in 
i,j € {0,p, r}. We also define •y ab = 7^ and 7 = det7 a b. 

The DBI action for our Dp/Dq system is given by equation ( |2.1| ), and can be ex- 
pressed as 2 

S DB i = -N f T Dq J d^xdrdOne-tJ-^gfaW, (2.8) 

where the coefficients g a b are those of equation ( |2.5| ) and 

W(r) = det 7lJ = g„ (g m g u + (2ita') 2 E 2 ^ + (2ita'f ( 9li A' 2 + g 00 A'*) . 

There are two conserved quantities, n q and j p , related to the charge density and the 
electric current respectively and defined by 

SSpBI _ nq_ SSpBi _ jp 

SA' Q n 3A' p Q n 

where Q n is the volume of the unit n-sphere. The gauge field components can then 
be solved for as follows 



A' - - Uq 



g rr goo (9oo9u + E p 2 (2rca') 2 ) 



gu \ Af 2 e-^g 00 g^ n + (2naf (g on 2 + 9ii j 2 ) 



(2.9a) 



n q g o 

where Af = NfTr> q fl n {27ia') 2 . As goo < (provided we are outside any possible 
horizons in the geometry), there is a radius r*(E p ) such that googu{ r *) + E p 2 {2Tia') 2 = 
0. This characterizes a hypersurface of constant r = r*. At this point the legendre 
transformed action vanishes on-shell. We will refer to it hereafter as the "singular 
shell". Calling the radius of the horizon, one can easily see that r*(E p ) > r^, and 
r*(0) = rh- The Legendre transform of the DBI action is Sbbi ~ — / d p+l xdrCoBi 
with 



r -<P n 7 10 ,\2m\ p- 1 ~ h j_ ( 27Ta ') 2 (9oon 2 q + guj 2 ) 

Cdbi = ~e *y/-{googu + (2vra') Efygfa g rr gnn\ 1 + K r2 p -24> n n v n n ■ 

v V J 9oo9u9un 

(2.10) 

Furthermore, as the quantity googu + E p 2 (2na') 2 changes sign at r = r*, in order for 
( p,10|) to remain real, we impose that the second square root in it also changes sign 



2 we will consider only fluctuations which are independent of the Wess-Zumino term. 
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at r*. This leads to a relation 3 between j p and E p : 



Jp = E p JX^gfr 2 g^ + nmva'fg 



-2 



(2.11) 



and defining the conductivity as usual through Ohm's law, we get a = j p / E p . This 
is the conductivity first found in [1] and the result is calculated purely in terms of 
background quantities. Although compact, the notation should not hide the impor- 
tant fact that a is a non-linear function of E p (see ( |3.5| ) below). This dependence 
arises from the evaluation of induced metric components at the singular shell r*. By 
studying linearized fluctuations on top of these background embeddings we will be 
able to extract information about properties of the gauge theory, and also calculate 
the conductivity using the Kubo formula. Agreement with the conductivity found 
in the macroscopic Ohm's law calculation (eq. [2.11|) will be a good check of our 
methods (see section ^ for details). 



3. Specializing to the D3/D7 system 

In order to picture what is happening in the general Dp/ Dq intersection we discuss 
here the possible embeddings in the case of the D3/D7 solution. In future sections 
we will discuss fluctuations about these solutions so we now introduce the coordinate 
system which will be made use of throughout. We find it convenient at this stage 
to use the radial variable, u, defined as u = (r/j/r) 2 . The pullback metric of the 
-D7-brane onto the D3-brane background will then be given by 

ds 2 vr 2 T 2 t . 2 1NJ2 (l-^ 2 + 4w 2 (l-w 2 )^ 2 ) , 2 M „ 

(3.1) 

where u = Uf, = corresponds to the UV of the gauge theory (the boundary of AdS) 
and u = Uh = 1 is the black hole horizon. The background field if)(u) corresponds to 
the classical embedding of the D7-brane, and the gauge field, including transverse 
perturbations, is given by 

A = A (u)dt + (E z t + A z {u)) dz + A±(t, z, u)dx ± , (3.2) 

with x± = {x,y}. Here A (u) corresponds to the non-zero baryon density, A z (u) 
to the non-zero current in the ^-direction, and E z is the finite electric field in the z 



3 it is interesting to note that in [23-25], a remarkably similar structure to the singular shell was 
observed in the trailing string solution, related to the local speed of light. Clearly the situations of 
a single moving quark and the quark current studied here are extremely similar. One can show that 
the infinite mass limit of the present current coincides with the current of a single quark moving 
through the plasma. We would like to thank C. Herzog and J. Casalderrey for pointing out this 
interesting observation. 
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direction. In this paper we will deal only with the transverse fluctuations, captured 
entirely in A±(t,z,u) and defer the full investigation of the longitudinal and scalar 
modes for future study. The transverse sector decouples completely as we consider 
momentum in the direction parallel to the background electric field (see section |] for 
further details). It is convenient to define the dimensionless electric field e as 

e = ^KE z , (3.3) 
(vrTi? 2 ) 2 V ; 

and constants of motion n q and J z , related to the charge density and current given 
in the previous section by 



n q \ __ 4(27TQt / ) / n q 



hj N c NfT 2 (irTR 2 ) Vj 



(3.4) 



However there is only one value of the current which renders a well defined action 
inside the singular shell. The position of the singular shell in this case is given by 
u+ — 1/ yl + e 2 , and the current is found by demanding reality of the action across 
the singular shell, leading to 



n 2 



h 



VTT^(i-^) 3 + —^, (3.5) 

where ij)* = ip(u*). From this we can read off the electrical conductivity 

j z N c N f T~u N c N f T , , 

a = dIL = c 1 J _i = 6 I a . 3.6 

E z 4?r e 4?r v 1 

Using equation (|3.5|) we can express the background gauge fields (Ao(u), A z {u)) 
entirely in terms of n q , e and the probe brane embedding profile ip{u) 



A'(u 



2 {7rTR 2 ) 2 n 2 q (1 - u 2 (1 + e 2 )) (l - ^ 2 + 4u 2 (1 - u 2 ) f 2 ) 



A(2na') 2 (1 - ^j 2 ) ((1 - u 2 ) (1 - ^ 2 ) 3 - (j 2 - ^r) - u^n 2 , 
n z (1 — u i y 

From here, Aq(u) can be obtained by integrating from the horizon u = 1 with 
boundary condition A (l) = 0. The UV series solutions to these gauge fields are 
given by 

TR 2 

A = n - -^-j-hgii + ... , (3.8a) 
A z = r ^ r ~ Jz u + ..., (3.8b) 
H being the chemical potential canonically conjugate to the baryon density. 
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The background field solution ip(u) corresponds to the classical embedding of the 
-D7-brane. The range of this field is < ip(u) < 1 and it is completely determined by 
its value either on the singular shell, or, for "Minkowski embeddings", the value of 
u for which this field goes to 1. The asymptotic behavior of the embedding near the 
UV gives the usual source (mass) and vev for the quark bilinear operator, leading to 
the relation that 

^(u) u ^ * ^ + . (3.9) 

The quark mass and condensate being given by 

1 



M n = ^\/9Y M NcTm g , (3.10a) 



M) = -y g * M N c N f N c T% . (3.10b) 

Although this is really the source and vev of a supersymmetric operator [26], in the 
present case where supersymmetry is broken, both by the electric field and by the 
finite temperature, we will simply refer to the operator in question as the quark 
bilinear. 

In what follows we shall plot the classical solutions to the embedding equation in a 
cartesian coordinate system in which it is more intuitive to picture the flows of the 
-D7-brane. The change of variables between the original coordinates and those we 
will use to plot the solutions is 



(p,l) = (VR?iW) 



1 - u 2 + 1 



u 



Notice that ip = cos 8 and the global factor defines an isotropic radial coordinate, 
w(u) = \/ p 2 + L 2 , as used in many papers (see for example [22]). 
For the D7 embedding in the finite temperature D3 geometry with finite electric 
field there are three families of solutions illustrated in figure |l|. We discuss here the 
three types of behavior labelled "Minkowski solutions" , "Black hole solutions" and 
"Conical solutions". 

Minkowski solutions: These solutions do not intersect the singular shell and 
therefore do not meet the horizon. The current in this case is zero as there are no 
free charge carriers. For these solutions the energy supplied by the electric field is not 
sufficient to pair create quarks of such a high mass. They therefore only appear for 
sufficiently high values of m q /T. Fundamental matter remains bound in long-lived 
mesons. 

Minkowski solutions are labelled by the point of closest radial approach to the D3 
branes, u = u max , at which they reach 9 = or ip(u max ) = 1 with v>h > > u max . 
The second boundary condition needed to fix these solutions determines that there 
is no conical singularity at u max , corresponding to ip'{u max ) = — oo. 
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h q =0, e=3 



4 - 




Figure 1: Embeddings at zero baryon density showing the Minkowski embeddings, which 
do not hit the singular shell (in this case at p 2 + L 2 = 3 + y/W) conical embeddings given 
by short dashed lines which pass through the singular shell and hit the L-axis, and black 
hole solutions which pass smoothly from the singular shell to the horizon. 



Black hole solutions: These solutions pass smoothly through the singular shell 
and end on the black hole horizon without ever reaching i/j(u) = 1 for any value 
of u. The fact that the solutions pierce the singular shell indicates that there is 
a finite current. For these solutions, the quark mass is sufficiently small that the 
electric field can pair-create them. With these free charges a current can be set up. 
With an external electric field there is an unlimited supply of energy and one might 
have thought that the charge carriers would continue to accelerate. However, in the 
presence of the adjoint matter (with Nf/N c —>■ 0) the energy of the fundamental 
matter is lost into this background and an equilibrium situation with a fixed, finite 
current is set up [1] . At finite baryon density all solutions have to end on the horizon 
(figure 0). 

Conical solutions: These pass through the singular shell and appear to display 
a conical singularity at a radius, u cs < Uh, determined by ip(u cs ) = 1. We illustrate 
one of these solutions in figure |3| in the original u, ip variables. 

Solutions which pass through the singular shell are determined uniquely by specifying 
ipi, = ip(u*). All derivatives are fixed by this single boundary condition and the 
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Figure 2: Embeddings at finite baryon density and the same electric field as in the 
previous figure illustrating the absence of Minkowski solutions. Note that there are still 
conical solutions in this case. 




Figure 3: Example of a solution with two conical singularities. The left hand plot includes 
the whole range of the embedding from the UV (u = 0) to the horizon (u = 1). The second 
part focuses just on the region around the conical singularity, illustrated in the left hand 
plot by the boxed region. Here the first conical singularity is labeled as (u cs ,^ cs ) while the 
second is labeled as (u CS 2,ip C s2)- The L-axis is the line ip = 1 while the horizon is marked 
as the black box. The singular shell is marked in the second plot as u*. 

singular shell acts as an attractor for solutions 4 . This is reminiscent of the black hole 
4 there are four solutions to the expansion about the singular shell, two of which are complex 
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horizon itself in the zero electric field situation, with which one can fix the behavior 
of a solution in the UV uniquely with a single boundary value. 
The electric field lines caused by the presence of free charges need to end somewhere 
and, hence, the brane must intersect the horizon. For the conical embeddings, the 
same reasoning applies. Solving numerically we find that these solutions continue 
past the point of the conical singularity until they reach the horizon. Still, in (p, L) 
coordinates, at the points where they touch the L-axis one can observe that there 
is a jump in slope, hence the name conical singularity. After reaching this point, 
the value of ip starts to decrease again, and the brane eventually hits the horizon 
(possible after encountering another conical singularity as illustrated in figure |3[). By 
studying the tension of the D7-brane at the conical singularity we find that there is 
no other object which can be attached to the D7-brane matching it, meaning that the 
physical configuration is the one where the D7-branes continues through the conical 
singularity, eventually reaching the horizon. 

3.1 Singular shell solutions 

The solutions which pass through the singular shell exhibit subtle behavior which will 
influence much of the physics in what is to come. For this reason we will now discuss 
the details of the singular shell embeddings, parametrised by the three variables m q , 
e and n q . 

Especially interesting is the behavior in the (m q , e) plane. First we chose to fix the 
electric field and look at the solutions as a function of m q . 

In figure ^ we see the behavior for a range of masses. For the largest mass, m q = 5, 




p 



12 3 4 1234 1234 12 3 4 



Figure 4: Evolution of the embedding profiles with decreasing values of the mass m q = 
5, 4.5, 4 and 3 at fixed values of h q = 0.05 and electric field e = 5.8. After developing 
conical singularities the profile spreads suddenly for lower masses. 

we see that there is a narrow throat, which leads from a relatively straight D7- 
embedding, sharply down into the singular shell and subsequently the horizon. As 
and one of which gives non-physical, oscillatory behavior. 
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we decrease the mass to m q = 4.5 we see that the behavior changes and we obtain 
a solution with a conical singularity. At this point the brane has a smoother profile 
and the sharp throat has disappeared. For still smaller masses than m q = 4 we 
return to the pure black hole embedding and the conical singularity disappears. The 
profile is much smoother and there is a "wide throat" to the solution, rather than 
the previous sharp spike. Decreasing m q still further widens this throat slowly and 
monotonically. 

It should be noted that, with decreasing m q , the conical solutions start with the 
conical singularity at the point where the black hole horizon meets the L-axis (u = 
Uh). The position u cs of the first conical singularity (for there may be multiple) 
moves away from the horizon and then returns to it as the quark mass is decreased 
further. 

We see from this example that over a small range of masses there is a very fast 
change in behavior, from a sharp spike solution, through a conical solution, to a 
smooth, wide throat solution. This behavior is generic, though the region over which 
this fast transition occurs is dependent on both m q and e. Indeed not all values of 
n q and e will exhibit conical solutions. 

Keeping m q fixed and varying e gives a similar behavior, exhibited in figure [5[ We 
label the region where the behavior changes quickly as a function of m q or e the 
"crossover region" . As stated above, this region may or may not contain a subregion 
of conical singular solutions, depending on the value of n q . 

L L L L 



5 h 5 h 5 h 5 h 




1234 1234 1234 12 3 4 



Figure 5: Evolution of the embedding profiles with the electric fields e = 0.8,3,5.8 and 
6 at fixed values of h q = 0.05 and quark mass m q = 4.5. We observe the developement of 
conical embeddings in the range e 6 (5.7,6). Increasing the value of e there is a sudden 
change in form, and the profiles spread significantly. 

As a physical motivation for such a crossover it seems that there are two competing 
mechanisms forcing the brane embedding to go in two different directions. The 
crossover is the point where one of of these effects becomes dominant over the other 
and there is a sudden change in the behavior. This is rather reminiscent of the 
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two effects seen in the expression for the conductivity, one of which becomes more 
important as e increases and one of which becomes more important as it decreases. 
These two terms are related to the the pair-creation and presence of free charges 
respectively. 

We can plot the crossover region in the (m q , e) plane in which conical solutions are 
present. This wedge is shown in figure |6] for varying values of n q corresponding to 
the increasingly dark colors. It is rather peculiar that the only effect of the baryon 
density is to shift the starting point within the same region and not to displace the 
region at all in some direction. 




Figure 6: The shaded region is where conical embeddings are found in the plane (m q ,e). 



The regions for different h„ are superimposed, beginning at larger m„ the larger h q 



is. 



The 



lowest region is for h q = 2 • 10 and this contains all other regions for larger values of n q . 
The orange region is for h q = 0.6 and the brown one corresponds to h q = 1.2. 



3.2 Chiral condensate 

Having detected a crossover, given by a sudden change in the behavior of probe 
embeddings, it is natural to expect that other physical quantities will be affected. In 
the present section we will examine the chiral condensate. 

We plot — c ~ (^/^/) given by equation ( |3.10| ) as a function of m q and e for n q = 0.012 
in figure [7[ We see that in the region where the crossover in behavior occurs at the 
singular shell the condensate also undergoes a faster change in behavior, as compared 
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Figure 7: Condensate as a function of electric field and mass for n q = 0.012. Curves of 
constant m q are plotted in blue, and those of constant e in green. We superimpose in black 
the curve where conical embeddings occur. Upon crossing it there is a mild but noticeable 
change in slope. 



with other regions in the (m q , e) plane. We have superimposed on the graph the curve 
e = 0.29 m 2 q which runs through the middle of the orange region in figure |6|. We see 
that this curve coincides roughly, as expected, with the region where the level curves 
change from concave to convex. We should note however that the change in behavior 
of the chiral condensate over the crossover region is much less pronounced than that 
of many other quantities as we will show. 



3.3 Charge susceptibility 

Another physically relevant quantity showing a fast change in behavior 
crossover described 
namically defined as 



, _ - across the 

crossover described in the previous sections is the charge susceptibility, thermody- 



dfx 



N c NfT 2 dn 



N c N f T 2 ~ 



(3.11) 



where /i is defined 5 in equation ( |3.8a| ). From eq. (|2.9a| ) we can read the dependence 
of \i on the charge density at fixed temperature and extract the value of the sus- 
ceptibility. In figure § we plot this quantity as a function of the quark mass and 
the electric field for n q = 0.6. There we can see that in the region where we have 



5 where we have defined An 



4a' v " 



fl q U + 



from equation (3.8a 
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narrow throats in the embedding profiles the susceptibility varies very slowly. On 
the other hand, in the region with wide-throat profiles we have a more pronounced 
slope. These two regions are separated by the line where conical embeddings may 
occur. We will find this behavior again, when considering the conductivity and the 
photoproduction. 




Figure 8: Dimensionles charge susceptibility H as a function of the quark mass and 
electric field for n q = 0.6. Green, blue and black lines have the same meaning as in the 
figure for the chiral condensate. 



3.4 Conductivity 

In [1], by analyzing eq. ( |3.6|) it was conjectured that the two terms present under the 
square root have physically clear and distinct origins. The second one is the clearest, 
as it is directly related to the presence of charge carriers proportional to the net 
baryon charge ii q . Notice the presence of the electric field in the denominator, which 
suppresses this contribution for larger values of e. The first one should correspond 
to the presence of charge carriers that come from pair production from the vacuum. 
It is enhanced by a stronger electric field, and, assuming a direct relation between ip* 
and the mass, this term should vanish for high enough masses. Roughly speaking, 
this behavior is seen in the set of graphs shown in figure || 

What is less obvious is the sharpness of the transition to the region of high conduc- 
tivity dominated by the first term in (|5.5|) . This is the slope which is visible for high 
values of e and low values of m q . Clearly the abruptness in this transition comes 
from the hidden dependence of ip± on e and m q , and the sudden change in profiles 
that was investigated in the previous section is clearly behind this strong sensitivity. 
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Figure 9: Dimensionless conductivity, a, as a function of e and m q for two different 
values of h q = 0.05 and 1.2. Varying m q at constant e shows a stronger change in slope 
upon crossing the black line than varying e. 



3.5 Charge diffusion 

A formal study of the charge diffusion for this system would involve solving for 
the longitudinal vector channel, that presents a pole in its correlator which, in the 
hydrodynamic limit, gives information about this transport coefficient. In this paper 
we do not follow this procedure but exploit the Einstein relation, which expresses 
the diffusion as the quotient of the conductivity by the susceptibility 

„ a 1 a D , . 

D = - = ^ = . 3.12 

H 2nTE 2vrT v ; 

With our numerical methods we are able to recover the vanishing electric field results 

of [27], and in figure [Tl] we show how these get modified in the presence of a finite 

electric field for h q = 0.6. 

In this figure we can see an interesting behavior. Keeping the electric field fixed and 
studying the diffusion constant when the quark mass increases we first see that it 
decreases very slowly until one reaches the crossover region. Once this happens there 
is first a rapid decrease and, once we pass this region and enter the parameter space 
where we have narrow-throat embedding profiles, we find that the diffusion constant 
increases with the quark mass very fast. This rapidly changing behavior with the 
mass is correlated with the vanishing of the susceptilibity in the same limit. 

4. Fluctuations 

In order to compute current- current correlators we must consider perturbations of 
the gauge field A^. The components Aq and A p are scalars under the little group 
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0(p—l), and will mix with any other scalar modes present in the system, for example 
the profile fluctuation 5ip or fluctuations of the dilaton 5(f). This channel is fairly 
complicated to analyse and for the D3/D7 case it has been studied in the absence 
of electric field in [19]. On the other hand transverse fluctuations A± = A. t with 
i = 1, ■ • • ,p — 1 will transform as vectors, and hence decouple completely from all 
other excitation modes and have been studied, again in the absence of an electric 
field, in [28]. In this paper we will be concerned exclusively with these transverse 
fluctuations, leaving the coupled sector for future study. 

In a rotationally invariant theory, the retarded Green's function G^ u involving per- 
turbations of the form exp(— ik^x 11 ) is generically split into two different polarization 
tensors [29] 

G>(*) = Pt,(k)IL L (k) + Pl(k)IlHk) , (4.1) 

where the projectors are defined as P q = P „ = 0, P^ n = 5 mn — k m k n /k 2 , pj}„ = 
Tijiv — k^kv/k 2 . In the present setup we expect this structure to remain unaltered 
as long as the background electric field, E p , points in the same direction, x p , as the 
propagation of the perturbation and, hence leaves the 0(p — 1) invariance untouched. 
For the transverse fluctuations, the components of G^ s for r, s — {1, • • • ,p — 1} are 
related directly to II 1 - by 

G*{k) = S^ik) . (4.2) 
The prescription for the Lorentzian correlator relates the polarization function n- L (/c) 
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to the boundary action for the transverse fluctuations 



S bou = ~ I d p+l x e-*v^7^7 rr A^± , (4-3) 



9 

Jr->r b 

with j\f defined as in equation ( |2.9|) . This prescription gives [30] 

A' A* 

n ± =ATe-* v ^ 7 ^7 rr 1 1 



\A\< 



(4.4) 



In section || we discuss further how to calculate this quantity for the specific case of 
the D3/D7 system. 

4.1 Transport phenomena 

The conductivity is defined macroscopically via Ohm's law as the ratio between the 
current density and the field that creates that flow of current. The calculation of the 
conductivity for a Dp/Dq intersection has been done following this definition in [1]. 
Here we are interested in a microscopic derivation using the Kubo relation: 

„ Imil- L (a; = q) 

= - hm — , (4.5) 

where D is the diffusion constant and H the susceptibility of the medium. In the next 
section we will show that this prescription leads to the same analytical expression 
obtained in [1], thereby proving the validity of the Einstein's relation a = DS. An 
independent check of this relation amounts to calculating the susceptibility 5 and 
the diffusion constant D separately. 

The susceptibility is an equilibrium quantity given by a thermodynamic definition 
S = dn q /dfi , where /i is the chemical potential and n q the charge density, as has 
been discussed in section I3TBI. On the other hand, the diffusion constant, D, can be 



obtained from a hydrodynamic pole in the correlator for the longitudinal channel 
of vector fluctuations, relating the frequency and the momentum as dictated by 
Fick's law: uj = —iDq 2 + 0(q 3 ) , although we leave this calculation for the future. 
Numerically one can see that in the case of vanishing electric field the product WE 
obtained in this way agrees with a, as explicitely shown in [27]. 

4.2 Gauge field fluctuations 

At this point we must study what boundary conditions we need to impose on the 
transverse fluctuations at any singular points in their equations of motion. The 
equation of motion for the transverse fluctuations after inserting the plane wave 
ansatz A±(t, x, r) — > e~^ UJX °~ qxP ^A±(r) reads 

d r [e'^V^ll 11 - (Y r A' ± - i (cu 7 {0r) - g 7 (pr) ) A ± )] - ie"^^^ (cu 7 (0r) - g 7 (pr) ) A' ± 
-e-^v/=7~ 7 ^ (^ 2 7 00 + gV P ) A± + 2 uj q e^v/=7~7^ 7 (0p U ± = . (4.6) 
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At finite temperature we need to make sure that the solution is regular at the 
horizon 6 . This can be achieved by performing a Frobenius expansion. The solution 
is regularised as A ± (r) = (r — r h ) v h Aj_ treg where the indices are found to be 

CO 1 ± 1 • u ( A 1\ 

V+ = i — — , 4.7 

'± 2 4vrT' v ' 

the (h) indicating that this is an index at the horizon. Note that this is markedly 
different from the case of zero electric field where the only possibilities were outgoing 
or incoming boundary conditions. Now there is the possibility of an additional damp- 
ing factor at the horizon. In this case the choice of outgoing or incoming boundary 
conditions is dependent on the relative sign between A' and A' p . In our choice of the 
background gauge fields (equation (|2.9| )) we have chosen the relative sign of these to 
give the appropriate incoming wave condition. 

In addition to the singular behavior at the horizon, there is also a singularity in the 
equations of motion for the fluctuations at the singular shell, r*. We can perform a 
similar Frobenius study in this region. In this case we find one trivial index, rj^ = 0, 
while the other, T]r*\ is a very complicated function of E p: u>, n q and ip* which, in 
the limit of vanishing electric field, is given by rjr(E p — > 0) = i-^j,- 
Finally, in order to choose our boundary conditions we want to recover the usual 
expression for the indices at the horizon in the limit of E p — > 0. Notice that in this 
limit, as we have already seen, r* — > r^, so 

(r — Th) = lim (r — nj (r — r+) = lim (r — nj = (r — ru) i7rT ■ 

(4.8) 

We have four possible combinations of the two values at the horizon and at the 
singular shell, combining to give the following values for r^( lim ). From table we can 



„(Hm) 


(h) 

v+ 


(h) 

?7_ 


W 

m 


1 - i 


1 A-kT 
1 AttT 



Table 1: The four possibilities for the index on the horizon, rf* ltm \ in the limit of vanishing 
electric field. 

see that to perform our study with a consistent E p — > limit we have to take the 
horizon index rj^ and at the singular shell 77Q * . 

As a further check that we chose the correct index at the singular shell we will show 
that our choice is the unique one that recovers the conductivity fl2.11|) via the Kubo 
formula 04.51). 



6 provided a non vanishing n q the branes can be black hole-type or conical-type. They touch the 
black hole horizon at some point, and the formula obtained will be shown to coincide with numerics 
also for conical-type embeddings. 
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4.2.1 Fluctuations in the hydrodynamic approximation 

Here we are interested in solving equation ( |4.6|) in the limit of vanishing momentum. 
It is a consistency check that the solutions are regular both at the horizon and at 
the singular shell. 

In the hydrodynamic limit we can expand the transverse fluctuation in a series in 
small lj and q. In order to perform this expansion we set (uj, q) — > \hyd{&, o) an d 
perform the expansion for small \h y d- The series for the fluctuation then reads 



Aj_ = A ( P + X hyd A ( l 1 + 0(X 



(i) 



hyd) ' 



(4.9) 



With the indices chosen in the last section we can write down the regular functions 
in the bulk order by order. At leading order in the hydrodynamic limit we calculate 
the equation of motion for the transverse gauge field fluctuation (|4.6|) and find 

0. (4.10) 



■*V=7T i± 7 TP ^f )/ ) 
This has a closed solution given by 



Af{r)=C + C l 



dr 



'7, 



' — ryy- 1 — 1 -^' 



(4.11) 



where the integral is taken from the horizon out to a radius r. This integral is 
divergent at the singular shell and so the only regular solution is the the constant 
one A { °\r) = C . 

Calculating the equation for the fluctuation at next to leading order we obtain the 
following equation 



'*V=7T il 7 ,T ^ )/ ) -^(e-^ v /3 77^7 (0r) C ) / 
which can be integrated twice to give 

A^(r) = 







— dr + C 2 . 



n, 



(4.12) 



(4.13) 



This expression is completely regular at the horizon but needs to be regularized at 
the singular shell. With this regularity condition we obtain 

jpCo 



Ci 



-zu 



(4.14) 



We now calculate the conductivity from the Kubo relation (|4.5| ), which relates the 
transport coefficient to the transverse correlation function. Using this prescription 
we obtain: 



a 



-M lim Im 



G\ + iue 



-77 7 



(0r) Cn 



ujC 



T'bou 



(27ra') 2 E p 
9oo9u 



1 

V'boh 



(4.15) 
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The second term vanishes when evaluated at the boundary and we recover the result 
of [1]. This is a good clue that we have chosen the correct index at the singular shell. 
Had we chosen the index rjr we would have found a different answer (incompatible 
with the result of the macroscopic conductivity calculation), which in the limit of 
vanishing electric field reads (Je p =o = —^ m E p -*o j P /E p , exhibiting a change of sign 
that can be traced back to the time reversal that also appears in table [TJ. With this 
positive result in hand we can now go on to calculate the spectral function away 
from the hydrodynamic limit which will ultimately give us the photoproduction of 
the strongly coupled plasma in the presence of an external electric field. 



5. Fluctuations in the D3/D7 system 

Here we specialise once again to the D3/D7 system and turn to numerical calculations 
of physically relevant quantities. Having shown that our choice of indices is the 
unique one to recover the conductivity while also recovering the correct E p — > limit 
(the usual incoming wave boundary conditions) we are now able to study the spectral 
function for mesons dissociated both by the electric field and by finite temperature 
effects. 

Using the radial coordinate introduced in section £| the gauge field can be expanded 
in the UV (small u region) in terms of the two independent solutions 

A± = BiA {1) (it) + B 2 A {2) (u) , (5.1) 

where 

oo 

A® {u) = 1 + J2 + h A {2) log u , (5.2a) 

i=2 

oo 

A®(u) = u + ^afV. (5.2b) 

All the coefficients a[ and h are determined by recursion relations 7 , where B m are 
the two connection coefficients. Applying the lorentzian AdS / CFT prescription (eq. 
we find for the two-point function 

N C N,T 2 Al'(u) _ N C N,T 2 ( B 2 



n± = -^r- H -tM = 2 U + h (1 + los( " ^ 0)) J ' (5 ' 3) 

leading to an expression for the transverse part of the spectral function 

X x = -2 Im n x = —N c NfT 2 Im = N c N f T 2 x^ ■ (5.4) 

B\ 



7 for example h = q 2 — to 2 with to = ui/(2ttT) and q = q/ (2ttT). 
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This of course will be a function of e, h q and m q leading to a large parameter space. 
We split this analysis into two separate areas: The first is confined to solutions with 
zero spatial momentum which will allow us to see the appearance of quasiparticle 
like peaks in certain limits of the parameter space, and also to compute the con- 
ductivity numerically as the slope a = 1/2 lim^o X ± ( UJ : 0)/u. The second will be 
solutions with lightlike spatial momentum in which we will be able to calculate the 
photoproduction rate. 

To perform the study we will numerically integrate equation (|4.6| ) from the singular 
shell to the boundary. As mentioned in the previous section, we will use a regularising 
index t]^ = 0, so that the function A±(u) is regular at the singular shell. Once the 
numerical solution is obtained we read the asymptotic behavior near the boundary 
and from that data we extract the coefficients B m , which will depend on all the 
parameters considered. 



5.1 Transverse fluctuations at zero momentum 

Setting q = allows us to plot the spectral function as a function of vo and examine 
the evolution of it across the parameter space spanned by (e, n q ,m q ). A fairly im- 
portant change in shape is observed upon crossing the conical embedding transition 
region in the (m q , e) plane (see fig. |6|). In figure [11] we show an example sequence of 



spectral functions for a range of values of e and with n q = 0.012 and m q = 3.5. Over 
a very small range, 3.4 < e < 3.7, the magnitude of the spectral function drops by 
around two orders of magnitude and the sharp peaks disappear completely, leaving 
a smooth profile with very wide oscillations which can no longer be interpreted as 
quasiparticles. This rapid crossover in behavior is a relic of the very fast change in 
behavior of the embedding in a small region of the (m q , e) parameter space. 
A physical observable related to the spectral function at zero momentum is the con- 
ductivity. A perfect agreement between our numerical computation of this quantity 
and that performed analytically, eq. (|4.15|) , is found 



5.2 Lightlike momentum and photoproduction 

Here we generalise the previous calculation to the case of lightlike momenta. This 
allows us to study photoproduction and discover how an electric field affects the 
brightness of the plasma [29, 31]. 

The emission rate for real photons is controlled by the spectral function evaluated at 
lightlike momenta u — |k|. Given that on the light cone IT|| = (otherwise we would 
have a divergence in the retarded Green's function), we can express the emission rate 
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e=2.5 



e=3.4 




w 



w 



Figure 11: The difference, Ax , between the spectral function \ (to) and the zero tem- 
perature limit with increasing values of the electric field e at fixed value of the baryon density 
(h q = 0.012J and mass (m q = 3.5J. After crossing the line of conical embeddings the peaks 
broaden very quickly into wide oscillations which cannot be interpreted as quasiparticles. 



entirely in terms of II_l as follows 



dk 3 e| M Imn ± (^ 



,V) 3 2|k| e"/ T -l 



(5.6) 



where e EM is the electromagnetic coupling constant. Since our system lacks isotropy, 
strictly speaking, the above expression only applies to photons with = u(l, 0, 0, 1) 
which propagate parallel to the electric field E = (0, 0, E z ) 8 . In figure ^ we see a 
series of curves that represent the value of dT(to)/dto normalized by (2%T)e EM /n 2 . 
The Boltzman factor means that the relevant structure can only be seen in the range 
of frequencies vo 6 (0,2). The curves are plotted for m q = 3 and increasing values 
of the electric field e = 2.3, 2.5, 2.7 and 3. A strong enhancement is observed upon 
traversing the crossover region. 

The rates are similar and exhibit a peak at to ~ 0.25. We could choose this as a 
representative value of the overall magnitude of the photoproduction, and see how it 
changes as we move across the (m q , e) space. Alternatively we can integrate ( |5.6|) in 
uo G (0, oo), which gives the total number density of thermal photons emitted parallel 



i we would like to thank David Mateos for pointing out this subtlety 
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Figure 12: Photoproduction density rate in the direction of the electric field, from low to 
high, at values of e = 2.3, 2.5, 2.7 and 3 with n q = 0.05. 



to the electric field (per unit momentum space solid angle). The result can be seen 



in figure 13. A remarkable similarity with the plots in figure is observed. Of course 




Figure 13: Integrated photoproduction (2-k) T / e EM as a function of(m q ,e 



a knowledge of the differential photoproduction rate for q not parallel to E would be 
extremely interesting. We leave this analysis for a further study. 
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6. Discussion and conclusions 

In this paper we have studied in detail the effect of an external electric field on the 
physical properties of a large N c finite temperature Yang Mills plasma with a single 
flavour in the quenched approximation. In addition to the strength of the electric 
field the free parameters of our model have been the baryon density and the quark 
mass and we have explored this three-dimensional parameter space thoroughly. The 
singular shell which is a direct consequence of a large enough ratio of the electric 
field strength to the quark mass appears to play a key role in the quantities of the 
gauge theory living on the AdS boundary. Although the singular shell does not have 
the geometrical properties to label it a horizon, in many ways its effects appear to be 
similar to those of a real black hole horizon. In particular the attractor-like behavior 
whereby the solution of the embedding is fixed uniquely by a single boundary condi- 
tion constrains the solutions greatly. Continuing the embedding inside the singular 
shell gives us solutions which may or may not have a conical singularity, depending 
on the values of e and m q for a particular value of n q . Indeed this conical singularity 
appears and then disappears as we traverse the (e, m q ) parameter space. The region 
in which this transition happens appears to play a key role as a crossover region in 
the behavior of many physical quantities. However, it may be interpreted that this 
region of conical singularities is only a symptom of what is happening at the singular 
shell. 

In [19] it was argued that the size of the effective horizon induced on a probe D7- 
brane greatly affected the spectral function, with a small induced horizon leading to 
a series of sharp peaks in the spectral function and long-lived quasiparticles, whereas 
a large effective horizon led to a smooth spectral function and the disappearance 
of the quasiparticle interpretation as the poles in the quasinormal mode spectrum 
disappeared into the complex frequency plane. This phenomenology is related to the 
fact that the incoming wave boundary conditions at the horizon allow the absorption 
of energy into the system more rapidly if the effective horizon is larger. Quasiparticles 
have shorter lifetimes if, in the holographic picture, their energy can be dissipated 
into a larger effective horizon. 

Although the interpretation of the singular shell as any sort of effective horizon is 
unclear, we do appear to have a very similar phenomenology in the work presented 
here. For a large induced singular shell area on the D7-brane (corresponding to a 
small value of ■0*) the spectral function is smooth, while for a small induced singular 
shell area, the spectral function is sharply peaked. The crossover region is then 
traced to a small region in parameter space where the area of the induced singular 
shell changes very quickly from a small value to a large value, allowing, for instance, 
for several orders of magnitude change in the photoproduction rate. 

One may ask if this is not due to the effective horizon itself, however it appears 
that this is not so. As we pass through the cross-over region, the size of the effective 



- 26 - 



horizon does not increase monotonically, whereas it seems that the widths of the 
peaks in the spectral function, for example, do. If the widths of peaks in the spectral 
function were determined by the size of the effective horizon as in [19] one should 
see the same monotonic behaviour. For this reason it appears that the important 
IR quantity is the area of the effective singular shell, where again the incoming wave 
boundary conditions which completely fix the UV behaviour, seem to have a similar 
effect to the horizon in the e = case. In addition, because of the attractor-like 
mechanism at the singular shell, effectively the behavior outside the singular shell is 
completely shielded from what is happening inside. Any additional force that one 
applied to the brane inside the singular shell will not be felt outside. It is however 
true that the changes at the singular shell are reflected both in the UV physics 
and the behavior of the brane within the singular shell. For this reason the line of 
conical embeddings simply appears to be a good marker of the region where we have 
discovered such interesting effects. 

It should be noted that the calculation for the susceptibility, eq ( ft.llj ), does appear 
to rely on the behavior inside the singular shell. This is because we have imposed the 
vanishing of the time component of the gauge field on the horizon. A check that the 
susceptibility does not in fact depend on the region inside the singular shell would 
be to calculate it via the Einstein relation for which we would need to calculate the 
diffusion constant from the longitudinal correlator. We leave this investigation for 
the future. 

There may be hints to the physics behind the crossover region coming from the 
analytic expression from the conductivity. Indeed we can see that there are two 
terms in this expression, one of which is clearly related to the size of the induced 
singular shell area and increases with increasing electric field strength, while the 
other decreases with the electric field strength. As you increase e or 1 — ipl the 
dominant term in the conductivity will change from one to the other. Because m q 
and not -0* is the more physical parameter in our gauge theory the contribution of 
the first term is less clear when we plot physical quantities as a function of m q . We 
conjecture that the two competing effects which are very clear in the expression for 
the conductivity (namely pair-production and free charges respectively) may explain 
the rapid change in behavior across the crossover region. It would certainly be 
interesting to investigate the interplay between the two competing effects in more 
depth. 

It would be extremely exciting if one could test these predictions of a fast crossover 
region in the laboratory. However, it appears that the current technical limitations in 
heavy-ion experiments do not allow for the investigation of the quark gluon plasma 
in the presence of external electric fields. What may be more practical would be 
a condensed matter system modeled holographically, and demonstrating a similar 
phenomenology to the one studied in this paper. An obvious example would be the 
defect D3/D5 system. 
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